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ABSTRACT \ 


The accuracy of a magnetic field experiment on a space- 
craft is sometimes limited by the spacecraft field and not 
by the accuracy of the magnetometer itself. It has recently 
been suggested by Ness and collaborators that N - 2 magneto- 
meters on a radial boom may be employed to determine the 
first N-l multipole contributions to the spacecraft field 
and thereby improve the total accuracy of a magnetic field 
experiment. In this study the total error for systems of 
one to four magnetometers is investigated. The optimal 
magnetometer locations, for which the total error is a 
minimum, are found for given boom-length, instrument errors 
and magnetic field models characteristic for spacecraft with 
only a restricted or ineffective magnetic cleanliness program 
like Mariner -Venus-Mercury 1973. It is found that the error 
contribution by the magnetometer inaccuracy is increased as the 
number of magnetometers is increased whereas the spacecraft 
field uncertainty is diminished by an appreciably larger 
amount. Therefore, for the field models chosen and a boom-length 
less than 10m from the spacecraft center there is always a 
gain in accuracy compared with one magnetometer. For example, 
for a boom-length of 6m, instrument errors of O.ly and a rather 
conservative model describing MYM73 the total error bounds 
are 2.1y, 0.856y, 0.822y and 0.807y for one to four magnetometers 
respectively. A less conservative model yields 2.1y, 0.486y 
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0.379y and 0.355y, respectively. Thus the results of this 
investigation support the conclusion that the multimagneto- 
meter method is especially advantageous for observing 
geoastrophysical fields which are of the same order or 
smaller than the spacecraft fields at the end of a boom 
of technically reasonable length. 
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Introduction 


The accuracy of a magnetic field experiment on a space- 
craft is sometimes limited by the spacecraft field and not 
by the accuracy of the magnetometer itself. This problem 
has been discussed recently in two papers by Ness and colla- 
borators [Ness, 1970; Ness et al. 1971]. A recently suggested 
method to partly overcome this difficulty is to use N - 2 
magnetometers on a radial boom which yields an estimate of 
the first N - 1 multipole contributions of the spacecraft 
field and the ambient field and thereby improve the accuracy 
of the experiment [Ness et al., 1971]. A dual triaxial 
magnetometer system is being planned to be flown on the 
Mariner-Venus-Mercury 1973 mission (MVM73) . 

It is the purpose of this paper to present an investiga- 
tion of the important problem of optimizing the position of 
the magnetometers along a boom of given length to yield a 
minimized total error. We shall restrict ourselves to at 
most four magnetometers, which seems to be a practical limit 
due to weight, power and financial considerations. In section 
1 the error analysis is presented, in section 2 the analysis 
is applied to some illustrative cases and the overall results 
are discussed in section 3. 
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1 . Error Analysis 


We assume that a radial magnetometer boom is given, which 
defines a line on which we choose the origin of our coordinate 
system O'. If all magnetic sources of significance are inside 
a sphere of radius a around the origin, the spacecraft field 
on the boom can be represented [Ness et al., 1971] exactly by 

) 2 + t 

for r - a, (1) 


CO 



1=1 


where i = 1,2,3 denote the x, y, z-components in a suitably 
chosen frame of reference. Note, that we have chosen the 
to be the multipole fields at r = r , the outer tip of 
the boom measured from the origin in contrast to the usage 
by Ness et al. [1971]. This choice makes the formalism 
more compact. As already noted by Ness et al. [1971] one can 
see from equation 1, that it is advantageous to place the N 
magnetometers on one boom so that each multipole contribution is 
characterised by only three vector components instead of 21 + 1 
quantities as in the case of arbitrary locations, where l 
denotes the multipole order. For N > 2 this has the consequence 
that with a given number of magnetometers higher multipole con- 
tributions can be removed from the total error. In addition 
one magnetometer boom is technically superior to several booms. 
The bars on B and B„ denote that these are the exact 

dill y 1 

physical quantities without error. At a distance r, the total 
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field vector which includes the homogeneous ambient field 
to be measured, is given then by 



( 2 ) 


A magnetometer at position r would then have the following 
reading 


B , . = B . + B 

obs,i zo,i am,i 



(3) 


where B . is the instrument error of the magnetometer reading 
zo,i 

due to zero-offset, noise and quantization uncertainty. If we 

assume N magnetometers for each component, on the boom at 

positions r. (k = 1,***N), we obtain N equations for each 
x, k 

value of i = 1,2,3. 


B 


obs, i, k 


B . . + B 

zo,x,k am,x 



(4a) 


Without loss of generality we can choose the r. , according 

i, k 

to the following inequalities: 

i = 1,2,3 (4b) 


r, = r - i 
1 i,l 


v > r 

1.2 1,3 


. > 

> r. - a 
x,N 


It should be noted that in our formulation the positions 
of the mono axial magnetometers for a given k could be different 
for the three components and indeed all 3N monoaxial magnetometers 
may have no two radial distances the same. For all cases of this 
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study r is taken to be r . We now try to approximate 
1,1 ± 

the spacecraft field by the first N-l multipoles. 

If the magnetometers were ideal detectors, i.e., without 
error, we would have the following N equations to solve for 


the ambient field for each i: 


B . . . = B 

obs , 1 , k am 


•* * 5 ■-&) 


1 + 2 


k = 1, . . .N 


The unbarred quantities B . and B „ . can be considered as 

am , i ti 

estimates of the ambient field and the first N-l multipole 

fields, respectively. This set of linear equations can be 

solved for B . . It is convenient to introduce the 
am, l 

r i 

quantities p . , = : 

1,15 r i,k 


am,i det(p i ) 


obs, i, 1 
* obs, i, 2 
J obs, i, 3 


obs, i, N 


3 

P i,2 

3 

P i,3 


4 N+l 

P i, 2 ’ ' ' ’ P i, 2 


det (p.) 

l 


P i,2 

3 

P i,3 


P i,N 


N+l 
‘ P i,2 


. P 
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Here p. stands symbolically for the set p. p. _...p. . 

1 ■ 1 j X 1 ^ Z 1 j IN 

The determinants can be expanded with respect to the first 

column 


N 

det = ^ ^ (-l) k + 1 sub-det (k; p^) 

k=l 


(6c) 


For sub-det (k; p^) a relatively simple expression can be 
found by using some general algebraic theorems: 


N 

o n-1 


sub-det (k; p.) = 

P . ' 1 1 (p i - P . ) > o 

l , n -a l , n x , m 

(6d) 

1 n=2 

m=l ’ . ’ 


n^k 

m^k 



The inequality is fulfilled if inequalities 4b are satisfied. 
For n = k or m = k a factor one has to be set in the products 


defining sub-det (k; p^) in equation 6d. 

B . is given by 
am, l 


^ k+1 sub ^ det (Pi^ k > 

B . = (-1) 1 3 — — v- B . . . 

am, x / v . det (p.) obs,i,k 

k=l 1 

N 

y ; ^ik^i^ B obs,i,k 
k=l 


(6d) 


where 6d also defines the matrix (3. r . We are interested here 

lk 

in the error in the estimate B . , i.e. , AB. = B . - B 

am, x i am,i am, i 

Using equation 4 and some general properties of determinants we 
obtain: 
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>i - E <”i> 


zo, i, k 


+ S P ik (p i^ * B RN, i, : 


where 


- 2 +1 
B li P i,k 


This represents the exact equation for the total error. The 

first part represents the error due to the inaccuracy of each 

single magnetometer and the second part the error due to the 

residual spacecraft fields . , determined by the multipoles 

RN, i,k 

N 

of orders 2 to infinity. An important point can be noted from 

equation 7. The error AB^ in the ambient field components is 

completely independent both of the ambient field B . and of 

am, x 

the first multipole fields B. . , with l = 1,...N - 1. 

A rigorous error analysis would use the statistical 
properties of B . . and B . , to derive first an expected 

ZO y 1 9 K XvJN } 1 ^ K 

value for AB.. Since the error statistics of B . , contain- 
i zo,i,k, 

ing noise errors, digitization errors and zero-offset errors, 
is relatively well known to be symmetric around B . = 0, we 

ZO J X y K 

obtain: 


J i > X] P ik (p i^ < B RN, i, k > ’ 
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where <AB^> could be referred to as systematic error. For 

T . > we could for example take the values from the last 

RN, x, k 

mapping of the spacecraft field before launch. After deriving 

<AB > the expected value of B . would be 
x am , 1 

<B .> = B . - <AB . > (9) 

am , 1 am , l x 

We would then define the statistical error as the value of 
3 

Cy ' (AB^ - <AB^>)^) 5 = 6B which with a very low probability 
i=l 

g 

(for example 10 ) will be exceeded statistically. Such an 

approach is not possible at the present time, since the statistics 
of unpredictable changes of spacecraft fields is generally only 
very poorly known. 

We use, therefore, a much simpler approach. First, we 
note that the equations for the single components in equation 
7 are not coupled. We can therefore restrict ourselves to 
one component and henceforth drop the subscript i denoting 
the components. The problem of optimization can then be stated 
as follows: Let all the magnetic sources be included in a 
sphere of radius a and let a boom of length r^ from the center 
of the sphere be given. Furthermore let the uncertainty of a 
magnetometer reading |b | - B > 0 be given as well as a 

ZO , K K 

worst case magnetic field distribution along the boom. Then 
the question is: For what distances of the magnetometers from 
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the origin r (k = 1,2,...N) is the error in a given 

K 

component least? We shall consider the cases N = 1,2, 3, 4. 
From equation 7 we obtain 


! N 

y'' sub-det (p; k) |b , I 

zo,k' 

k=l 

N 

+ 1^ (-l) k+1 sub-det (p; k) B RN J 
k=l 

I N 

Z SUb ‘ det <P;k) ' B M,k 

k=l 

N 

+ (-l) k+1 sub-det (p ; k) B RN fc | 

k=l 


( 10 ) 


( 11 ) 


AB can be considered as an upper bound on the error AB. 
max 

2 . R esults of Optimization of Upper Bound on Error 

AB for Various Parameters and Field Models 
max 

In this section we apply equation 11 for several 

spacecraft field models and geometrical configurations. Before 

doing this two points may be noted. Because of the remark 

after equation 7 the quantities B giving the spacecraft 

RN j K 

field at r = r k without the first N-l multipole contributions 
could be replaced by the complete magnetic field components 
B i, at r = r T> Secondly the following "similarity laws" hold: 

SC « K K 
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r 


AB 

max 




B 


N; 


M, 2 


B 


M, 1 


B B 
M, N sc 

Vi 


T<) 


( 12 ) 


Thus, if the accuracy of all magnetometers is improved by 

a factor of q and the value of the spacecraft field is 

decreased by a factor q everywhere, the optimal magnetometer 

positions do not change, and AB is reduced by a factor 

’ max 

of q,etc. 

We shall present results for the following choices 
of our parameters: 

a = 1 m 

4m - r - 10 m 

B . = 0.1 v for all k 
M, k 

The choice a = lm for the sphere around the sources 
corresponds approximately to the main body of a Mariner-type 
spacecraft (without appendages) or an OGO-type spacecraft, 
whereas application to the smaller IMP- and Pioneer-type 
spacecraft would be possible using the similarity results. 

The origin O' will then be close to the geometrical spacecraft 
center . 

The value of 0.1 y for the absolute accuracy of a single 
magnetometer is achievable by present day flux-gate magneto- 
meters, with flipper devices for example, for a bandwidth of 
1 Hz and a digitization window of 0.02 y. 

The following field models are used: 


11 



B SC (r) - B dipole ( F )3 2] < 13 > 

1=0 

with the values of B,. , and f given in Table 1. 

dipole 

All magnetic models contain multipole contributions up 
to the 32-pole. The models of odd number have equal multi- 
pole contributions at r = a [see for example Ness et al, 1971], 
whereas the even numbers describe models with equal multipole 
fields at r = 2a, an even more conservative estimate of the 
nature of the spacecraft field. The success of the multi- 
magnetometer method depends on the relative size of the field 
components B RN ± compared with the total spacecraft field 
components B gc ± which determine the spacecraft field part 
in the total error for one magnetometer. For a given dipole 
field, B ( j i p 0 ^ e J the case f = 2 can be considered as the worst 
case in the sense that the magnitude of the multipole fields 
of order l - 2 is appreciably above the average fields for 
a uniform and isotropic distribution of dipoles in the volume 
of the spacecraft. This is shown in the paper by Neubauer 
and Schatten [1972], Specifically their analysis indicates 
that for an isotropic and uniform distribution of dipoles in 
a sphere of radius a and for the components transverse to the 
boom, the rms-dipole field exceeds the rms-quadrupole, rms- 
octupole and the rms-16-pole fields above distances of 1.90a, 
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1.67a and 1.54a, respectively. For the usually larger radial 
fields the corresponding numbers are 1.64a, 1.51a, 1.42a, 
respectively. Since in addition large portions of the sphere 
of radius a approximating the MVM73-spacecraft , which is of 
special interest here, are empty, the transition distances 
would be even closer to the spacecraft. In conclusion the 
field models 2, 4,6,8 are rather conservative in their assumptions 
concerning the smallness of the multipole fields. 

The models 3 and 4 give equal fields at r^ = 6m and 
thereby illustrate the influence of f = 1 and f = 2 on the 
multimagnetometer errors for a boom length of 6m and a given 
error of 2.1y for one magnetometer. 

Table 1 also gives the magnetic field of models 1 to 8 at 
r = 12 ft. These values are spread around the estimate of 12. 5y 
at 12 ft (contained in the NASA Proposal Briefing Material 
and assumed here to represent a magnetic field component) for 
Mariner-Venus-Mercury 1973. Note that models 3 and 4 are close 
to this estimate. The models therefore describe magnetic 
fields typical for spacecraft with a negligible or ineffective 
magnetic cleanliness program. 

The results for the optimal value of the error bound 

AB for models 1 to 8 are shown in figure la,b,c,d as a 
in 3.x 

function of r^ the boom length. Table 2 contains the results 
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for models 3 and 4 for the optimal magnetometer locations 

and the contribution to the errors by the magnetometer 

inaccuracies and the factors for the calculation of the 

ambient field for r^ = 6m. It should be kept in mind 

that the maximum error AB (p) applies to one magnetic 

max 

field component. The error for one magnetometer has been 

calculated by using B__ (r ) + B^ . 

SC 1 M 

3. Discussion of Results and Conclusion s 
We have solved the problem of minimizing the upper error 
bound for one magnetic field component for various field 
models and boom-lengths. The minimum is found as a trade- 
off between the error due to the inaccuracy of the magneto- 
meters and the error due to the unknown spacecraft field. 

The first error is decreased by increasing the distances 
between the magnetometers whereas the latter error is decreased 
by moving the magnetometers away from the high field regions 
close to the spacecraft or in other words closer to each 
other towards the outer tip of the boom. 

We first restrict ourselves to the models 1-6, in which 
all multipole fields have the same sign. It is seem from 
figures la,b,c, that AB max can be reduced quite appreciably 
by using several magnetometers. The biggest gain is achieved 
by going from one to two magnetometers. In the cases of the 

field of 2y at r 1 = 6m the upper bound, AB = 2.1y for one 

l max ' 
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magnetometer should be compared with ^ B max = 0.486y for 
model 3 and 0.856y for model 4 as shown in table 2. The 
results essentially confirm those of Ness et al. [1971 ]. 

Since for smaller spacecraft fields the inaccuracy of 
the magnetometers represents a larger fraction of AB , we 
would expect the gains to be smaller. This is confirmed 
by the curves for models 1,3,5 and models 2,4,6. Since with 
a given spacecraft field distribution, such as in models 
1 to 6, the individual multipole contributions fall off 
faster as their order increases and since the contribution 
caused by the inaccuracy of the magnetometers increases with 
the numb dr of magnetometers, N, the optimal number of magneto- 
meters decreases with increasing length of the boom r . 

For the field models 1-6 used in this study, which are 
representative for spacecraft with a restricted magnetic 
cleanliness programs like MVM73, the boom length at which 

AB (N = 1) = AB (N = 2) lies far beyond 10m even for 
max max 

the low field in models 1 and 2, as figure la shows. 

For magnetically clean spacecraft we find the transition 

distances characterised by AB (N = 1) = AB (N = 2) very 

max max 

close to the spacecraft so that reasonably short booms and 

one triaxial magnetometer can be used. The transitions from 

AB (N = 3) < AB (N = 2) to AB (N = 3) > (N = 2) 

max max max max 

and the transition from N = 4 being better than N = 3 to 
the inverse relations are shown in figures la,b,c. For the 
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low field model 1 and = 4m, N = 4 is the best choice 
among N = 1,2, 3, 4 while for r = 6m, N = 3 is best. For 
model 2 N = 2 is the best choice for all values of r from 
4 m to 10 m. 

As the spacecraft field increases, the value of r , at 

which the transition occurs, moves outward. Models 3,5,6 

show the ordering AB (N = 4) < AB (N = 3) < AB (N = 2) . 

max max max 

Whereas in the low field models and long boom cases not very 
much is gained in going from 2 magnetometers to 3 or 4 magneto- 
meters, in the strong field cases at smaller boom-lengths 
appreciable gains may be achieved in using 3 or 4 magnetometers. 
For field model 5 and = 4m a factor of 2 can be gained in 
going from 2 to 4 magnetometers. 

Another point should be noted in comparing different values 

of N. Even if AB (N = 4) ^ AB (N = 2) the greater 

max max 

effort with 4 magnetometers may be worthwhile, because in 
AB max (N = 4) the spacecraft field contribution is much 
smaller than in AB max (N = 2) (see for example table 2) . The 
error is then determined much more by the relatively well 
known error properties of the magnetometer than by the 
generally poorly known characteristics of the spacecraft field 
and its changes. 

So far we have not considered models 7 and 8, which are 
characterized by the same absolute values of the multipole 
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fields in models 5 and 6 but with alternating signs. The 

results of these models demonstrate two logically connected 

points. Firstly, for given absolute values of the multipole 

fields, the worst case for AB is the case of equal signs 

of all multipole fields. Secondly, in using certain field 

models for optimizing the multimagnetometer system, care 

has to be taken that the optimum is not too strongly dependent 

on the exact field characteristics. For example, the low 

values of AB (N = 2) in models 7 and 8 are caused by 
max 

B rn = 0 at r = 1.0m and r = 2.0m, respectively. This zero 
is given, however, as a difference of large numbers and is 
therefore very sensitive to changes in the spacecraft field. 

As the calculations of Neubauer and Schatten [1972] show, 
beyond approximately 3a, where a is again the radius of the 
sphere enclosing the magnetic sources, the nondipolar field 
is dominated by the quadrupole part for most possible models. 
The problem just mentioned could therefore easily be resolved 
by imposing the further condition ~ 2.5a for the case 
N = 2 and less restrictive conditions for the innermost magneto 
meter in the cases N = 3 and N = 4. 

Finally we note some important items which have not 
been considered so far: 
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1. The accuracy of the magnetometer closest to the spacecraft 

is not critical in the overall error bound AB . For 

max 

example, for model 3 with 4 magnetometers table 2 shows 
that even an error bound of l.Oy at r^ would change 
ABmax by only 0.03y. In addition, the optimization could 
also be performed with varying B . An an example we 

M y K 

mention the final layout of the MVM73 dual magnetometer 

system. Since only at r^ a flipper mechanism is used 

we have approximately B„ .. = O.ly and B = 0.3y. The 

M, 1 M, 2 

actual magnetometer locations are r = 6.95m and r = 4.65m. 

-L / 

For the dual magnetometer system we obtain values of 

AB ma x of °-41y, 0.67y, 1.04y and 0.43y for models 2,4,6, 

8, respectively. The corresponding values of AB for one 

max 

magnetometer are 0.43y, 1.22y, 2.4y and 1.37y respectively. 
Thus for model 4 an improvement by a factor of 1.94 and 
for model 6 by a factor of 2.3 is achieved, although the 
assumption of f = 2 can be considered to be rather con- 
servative. These results do not change appreciably by 
including the fields from the solar panels. The actual 

magnetometer locations are optimal for Bj. „ = 232y 

dipole 

and f = 2. 

2. For real, finite size magnetometer-sensors, errors could 
also be due to gradients in the spacecraft field. The 


radial field component for a sensor of length L, which 
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averages 


over the field along its axis as a flux-gate 
does, would have a gradient error of 

2 1 + 2 

■ (2 + l) (3 + t) ,L., p / 1\ 

24 ( F B e ( 

r i ' t + 2 

due to the multipole field B e (— — ) . Equation 5 

can easily be changed to take it into account, apart 
from the fact that is is very small in most practical 
cases. 

3. Equation 5 implicitly contains the assumption that 

approximating the spacecraft field by the first N-l 

multipoles is the best choice for the representation 

of the spacecraft field. This choice, as well as the 

field model employed, is dictated by the generally rather 

poor knowledge of the properties of the spacecraft fields. 

If some additional information is given in a specific 

case, it should be used to improve the accuracy i.e. to 

diminish the error bounds AB along the lines sketched 

max 

in section 1 between equations 7 and 10. The work started 
in the paper by Neubauer and Schatten [1972] can be con- 
sidered as a beginning in this direction. 

In conclusion we have shown that under a wide variety 
of assumptions about the spacecraft field, the use of two to 
four magnetometers can appreciably increase the accuracy of a 
magnetic field experiment with a magnetometer boom of reasonable 
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length on a spacecraft with significant spacecraft field 
contributions . 
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TABLE 1 SPACECRAFT MAGNETIC FIELD MODELS 


Model 

^dipole 

f 

B„ at 6m 
SC 

B gc at 12 ft 
3 . 65m 

1 

100y 

1 

0 . 56y 

2 . 81y 

2 

100y 

2 

0. 69y 

4 . 29y 

3 

360 . 0463y 

1 

2y 

lO.lly 

4 

289 . 1901y 

2 

2y 

12 . 40y 

5 

550y 

1 

3 . 06y 

15. 45y 

6 

550y 

2 

3 . 80y 

23. 59y 

7 

550y 

-1 

2 . 18y 

8 . 84y 

8 

550y 

-2 

1 . 92y 

7 . 62y 
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TABLE 2. CHARACTERISTICS OF MULT I MAGNETOMETER SYSTEM FOR r 1 = 6m (SEE 
ALSO TEXT), WHERE THE ERROR FOR ONE MAGNETOMETER IS LESS THAN 2.1^ 
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FIGURE CAPTIONS 

Fig. la-d. Upper bounds AB on magnetic field errors for 

max 

spacecraft field models 1-8 as a function of boom- 
length. N is the number of magnetometersi'employed. 
The model number refers to Table 1. 
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Figure la. 
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Figure lb. 
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